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The Tevatron has measured a discrepancy relative to the Standard Model prediction in the 
forward-backward asymmetry in top quark pair production. This asymmetry grows with the rapid- 
ity difference of the two top quarks. It also increases with the invariant mass of the tt pair, reaching, 
for high invariant masses, 3.4 standard deviations above the Next to Leading Order prediction for 
the charge asymmetry of QCD. However, perfect agreement between experiment and the Standard 
Model was found in both total and differential cross section of top quark pair production. As this 
result could be a sign of new physics we have parametrized this new physics in terms of a complete 
set of dimension six operators involving the top quark. We have then used a Markov Chain Monte 
Carlo approach in order to find the best set of parameters that fits the data, using all available data 
regarding top quark pair production at the Tevatron. We have found that just a very small number 
of operators are able to fit the data better than the Standard Model. 



I. INTRODUCTION 



The most recent measurement of the forward-backward asymmetry, A l p B , in top quark pair production 
at the Tevatron [l], 0] was performed by the CDF collaboration using a data sample with 5.3 fb _1 of 
integrated luminosity After background subtraction, the value of A*§, B in the center-of-mass (CM) 
frame of the top quarks is 

A% = 0.158 ±0.074 (1) 

which constitutes about two standard deviations above the Next-to-Leading-Order (NLO) Standard 
Model (SM) prediction [| 

Ap§ M = 0.058 ± 0.009 . (2) 

Despite the discrepancy in Ap B , the total tt production cross section is in good agreement with the SM 
prediction. In fact, with 4.6 fb _1 collected luminosity, the top quark pair production cross section Q 
yields the result 

C7^ easured = 7.70 ± 0.52 pb (3) 
for a top quark of mass 172.5 GeV, which is in good agreement with the theoretical prediction [|| 

off M (MCFM) = 7.45t°;g3 pb (4) 
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where MCFM stands for Monte Carlo for FeMtobarn processes 0] ■ Measurements of the ti differential 
cross section with the ti invariant mass {m t i), do j dm t i were also performed by the CDF collaboration 
With an integrated luminosity of 2.7 fb _1 CDF has tested the m t f spectrum for consistency with the 
SM. The results are presented in table HI They have concluded that there is no evidence of non-SM 
physics in m t i distributions. Hence, whatever new physics explains the forward-backward asymmetry in 
ti production, it has to comply with all other measurements that are in agreement with the SM. Finally, 
measurements of the asymmetry for two regions of the top-antitop rapidity difference (Ay) and for two 
regions of the invariant mass {m t i) were performed by the CDF collaboration in Q . The results are pre- 
sented in table [TT| together with the theoretical predictions . The asymmetry at high mass is 3.4 standard 
deviations above the NLO prediction for the charge asymmetry of QCD. Recently the electroweak con- 
tributions to the asymmetry were re-analysed flo| just to conclude that the observed mass-dependent 
forward-backward asymmetry still shows a 3c deviation in the high mass region. Inclusion of corrections 
beyond NLO does not change this picture as well [ll|. The separate results at high mass and large 
AF contain partially independent information on the asymmetry mechanism. Therefore, a total of 14 
observables were measured at the Tevatron. This set of experimental values will be used to investigate 
whether the complete set of effective dimension six operators is able to describe the possible new physics 
responsible for the observed discrepancies while retaining the measurements in agreement with the SM. 
Recently DO [l2j has measured Ap B with 5.4 fb _1 of collected luminosity. As discussed in 13], their 
analysis does not observe a significant rise of the folded detector level asymmetry with respect to m t f— 
and Ay. Until these results are unfolded they can not be directly compared to the CDF ones, even if 
at the detector level they appear to be consistent within errors. We did not use the results [l2j in our 
analysis. 



Bin a (CDF result) a (SM-NLO) 



(GeV) 


(pbarn) 


(pbarn) 


350-400 


3.115 ±0.559 


2.450 


400-450 


1.690 ± 0.269 


1.900 


450-500 


0.790 ±0.170 


1.150 


500-550 


0.495 ± 0.114 


0.600 


550-600 


0.285 ± 0.071 


0.400 


600-700 


0.230 ± 0.073 


0.310 


700-800 


0.080 ± 0.037 


0.100 


800-1400 


0.041 ±0.021 


0.036 



TABLE I: CDF measurements of do/dm t t (integrated in each bin). We bin- wise scale our SM result (at LO) 
to match the SM-NLO result to emulate a m t i dependent fc-factor for fitting. The SM-NLO values are extracted 
from the plot in |l4j ]. 



Observables CDF result SM prediction 



A%(\AY t \ < 1.0) 
Af B (\AY t \ > 1.0) 
A'pgirritt-) 
A t FB (m t t+) 



(0.026 ±0.118) 
(0.611 ±0.256) 
(-0.116 ±0.153) 
(0.475 ±0.114) 



(0.039 ± 0.006) 
(0.123 ±0.008) 
(0.040 ± 0.006) 
(0.088 ± 0.013) 



TABLE II: CDF measurements 0] and SM predictions for the Forward-Backward Asymmetry for two regions 
of AY t and for two regions of m t j. A t FB (m t f+) stands for J 4" B (m t j > 450 GeV), while A PB (m t i—) stands for 
Af B (m t t < 450 GeV). 

There have been several attempts to explain this discrepancy. The most popular collection of models 
among theorists when trying to account for the Tevatron results are the ones with new gauge bosons, and 
in particular, axigluons, W' and Z' bosons [l5|. Explanations in the framework of SuperSymmetric models 
were discussed in fl6| . Other possible justifications for the inconsistency between theory and experiment 
in the asymmetry, while leaving the cross section for ti production within measured uncertainties, include 
i-channel exchange of color sextet or triplet scalar particles [17j , s-channel coloured unparticle contribu- 
tions [18] or s-channel new colour octet vector bosons contributions [19}, light coloured particles from a 
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particular SU(5) GUT model 0, extra dimensions [Hj], SO(10) models [H[, SO(5) <g) U(l) gauge-Higgs 
unification models [23} , new heavy quarks [24|, diquark models [25j and models where SU(3) C QCD 
theory is extended to SU(N) C which is spontaneously broken at a scale just above the weak scale [2a ]. 

The search for resonances decaying into tt has also been carried out at the Tevatron [27[ (see also [28| ) 
with negative results. CDF has tested vector resonances with masses between 450 GeV and 1500 GeV 
with widths equal to 1.2 % of their mass. With 4.8 fb _1 of integrated luminosity they found no evidence 
of resonant production of tt candidate events. This result sustains the argument of integrating out 
new heavy fields and strengthens the idea of adopting a model independent approach in explaining the 
measured asymmetry at the Tevatron. An independent approach, with the recourse to higher dimension 
operators was already discussed in [2jj. In this work we propose to study the effect of dimension six 
flavour changing neutral current (FCNC) operators together with four fermion (4F) operators. In order 
to find the best set of parameters that fits the data we will use a Markov Chain Monte Carlo (MCMC) 
approach. 

The paper is organized as follows. The next section is devoted to describe the effective operator approach 
and the number of independent operators that will be used in the analysis. In section 3 we describe the 
MCMC method and we present the results for the complete set of operators. Finally a discussion on the 
results and the conclusions are presented in section 4. 



II. THE EFFECTIVE OPERATOR APPROACH 



The Standard Model of particle physics is the low energy limit of a more general theory which could 
manifest itself through a set of effective operators of dimensions higher than four. The effective operator 
formalism assumes that this more general theory would be visible at very high energies and at an energy 
scale lower than A, the set of higher order operators would be suppressed by powers of A. The Lagrangian 
of the new theory can be written as a series in A with operators obeying the gauge symmetries of the SM 
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where C SM is the SM lagrangian and and contain all the dimension five and six operators 
respectively. This formalism allow us to parametrize new physics, beyond that of the SM, in a model- 
independent manner. The term £^ 5 -* is eliminated by baryon and lepton number conservation. Thus, any 
new particle or interaction is hidden in the dimension six operators which are listed in [30l l3l| . 

We divide the dimension six operators in two groups, the four-fermion (4F) operators and the non-4F 
operators. The later can then be grouped according to the gauge boson present in the triple vertex. As we 
are discussing tt production, the non-4F operators contributing to the process have at least one top quark 
in the interaction. Operators with one top quark, a light up-quark and one gauge boson will be called 
FCNC operators. If the gauge boson is a gluon they are classified as strong FCNC operators [H, l33j ; 
otherwise they will be called electroweak FCNC operators [U, [35[ . 

When looking for new physics that would explain the tt asymmetry in the framework of the effective 
operator approach we start by looking at the dimension six non-FCNC operators. As the final state is 
it, the only possible new vertex is an anomalous gtt interaction. Its contribution to the process would 
originate from the diagrams presented in Fig. [TJ However, because the initial state is symmetric, these 
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FIG. 1: Feynman diagrams for tt production via anomalous gtt interaction, 
diagrams will only contribute to the cross section but not to the asymmetry. Therefore, any changes 
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produced by these operators would change only the cross section, but not the asymmetry where the 
discrepancy is. The next class of operators we discuss are the FCNC ones. In this case there are two sets 
of diagrams to consider: the ones initiated by gluons and the ones initiated by light quarks. For the same 




FIG. 2: Feynman diagrams with FCNC operators for tt production via gluon fusion. 



U, C 



U, C 




FIG. 3: FCNC Feynman diagrams for tt production via qq fusion. 

reasons discussed for Fig. [T] the contribution of the diagrams in Fig. [5] to the asymmetry is negligible. 
Therefore there are only contributions coming from the diagrams in Fig. |3l Regarding those diagrams 
(Fig. [3|), and taking into account that the contribution of the c-quark is much smaller than that of the 
u-quark, we discard all contributions that have a c quark in the initial state. Note that there are no 
s-channel contributions for the FCNC case because we have a top-antitop pair in the final state. Finally 
we will consider all 4F-fermion operators as shown in Fig. |4j Note that we will consider not only the 




FIG. 4: Four-fermion Feynman diagrams for tt production. 

interference term with the SM contribution, of order 1 /A 2 but also the modulus-square terms of order 
1/A . We will now discuss the minimum number of operators to be used in our analysis. 



A. Effective operators in the strong sector 



Following the notation of [30] , the operators contributing to the strong FCNC vertices can be written as 

0%4> = ti^^uUG^, (6) 

and 

0% = 4 A Q ^D v u R G^ , 0% = ql \ a ^D„qi L G^ , (7) 

where G" is the gluonic field tensor, u l R stands for a right-handed quark singlet and q l L represents the 
left-handed quark doublet. FCNC occurs because one of the indices is always equal to 3 while the other is 
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either 1 or 2, that is, there is always one (and one only) top-quark present in the operator; the remaining 
fcrmion field in the interaction is either a u or a c-quark. Throughout this section we assume that O lJ 
and O ji are independent operators and the hermitian conjugate of all the operators are included in the 
final Lagrangian. The operators in (Eft are related to the operators in (JjfJ through equations of motion 
that also involve 4F operators |30 |433 .. 36]. However, the 4F operators appearing in those equations have 
either one or three top-quarks [36|. Therefore, if those 4F operators can be discarded, operators in (J7J 
can be discarded as well. The operators presented in this section will give rise to the FCNC vertices of the 
form gtui (with Uj, — u , c) and the corresponding hermitian conjugate interaction with an independent 
coefficient. 



B. Effective operators in the electroweak sector 

There are also effective operators stemming from the electroweak sector that would give rise to new FCNC 
interactions involving the top quark [34], |35[ . We start by listing the chirality flipping operators which 
are the equivalent to the ones in the strong sector, the only difference being the gluonic tensor replaced 
by the U(l) and SU(2) field tensors. They can be written as 

0% = & ^ < ~4> B» v , &> W4) = <f L n ^ 4 4> W' v , (8) 

and 

0% = 4 7/1 A,«j, , 0% = ql ltl D v <f L , 0^ w = <? L rx 1„D V qi , (9) 

where B^ v and are the U(l)y and SU(2)l field tensors, respectively. There are also equations 
of motion in the electroweak sector that relate the operators in © with the ones in © and with 4F 
operators [36] . A similar analysis to the one performed for the strong sector regarding the contribution 
of the 4F operators leads to the conclusion that we can neglect the operators (0 in our analysis. 

Besides chirality-flipping operators there are chirality conserving FCNC operators. Their flavour con- 
serving counterparts are already present in the SM lagrangian at tree-level. In fact, the vertex itZ has 
two vector contributions of different magnitudes, one proportional to 7 M 7l and the other proportional 
to 7^ 7r. Hence the flavour conserving contribution would modify the Z boson neutral current. All the 
chirality conserving operators involve the Higgs doublet. As the Higgs field is electrically neutral, there 
are more effective operators which will only contribute to new Z FCNC interactions. This set of operators 
can be written as 

0% = itfDpfifa'fvti, (10) 
= * &D»<t>) (6 r & . = * ri D*4>) («1 T, ql) , (11) 

and 

0% u = (<t L D»v> n )D M t, 0% u = (D^qlu^D^. (12) 

Again, the use of the equations of motion allow us to discard the operators in (fl"2j) . In the electroweak 
sector, there are now 4F operators with one top and one anti-top. However, those 4F operators always 
have one b-quark in the interaction or, if not, are CKM suppressed making its contribution to the tt 
asymmetry negligible. Furthermore, as was shown in (36[, for all the operators in (fTU)) and (fTTj) . O 1 ^ 
and Qi % are not independent. This means that the number of independent operators in (|T0|) and ([TT]) 

is reduced to three (for each light flavour). Finally, for this particular study, we can group Ol^ and 

®4>q under the same Lorentz structure which further reduces the number of independent operators 
in (fTTj) to two for each light flavour. 

The above discussion leads us to the conclusion that the minimum number of operators needed to describe 
the asymmetry is 8 for each light flavour. 

The amplitudes were generated with CalcHEP (37j , the Feynman rules for the effective operators were 
derived with LanHEP [38| and the integration was performed using the CUBA library [3!| . 
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C. Four-fermion operators 



We now turn to the four-fermion operators. In order to make the analysis as clear as possible we will 
reduce the operators to a manageable number making use of all allowed reduction procedures, from 
equations of motion to Fierz identities. Again, because the largest contribution to ti production occurs 
in uu fusion, we will discard all non u-quarks contribution in our study. We end up with a total of 12 
operators in agreement with [ioj . that is, 12 operators for each light up-quark flavour and we do not 
consider operators with down-quarks in the initial state. This simplification allow us to find hints of the 
type of operators that can contribute to the asymmetry according to the its Lorentz structure. We write 
the four fermion lagrangian as 

A.B 

% E 1 N ab (u Al ,t A ){i Bl »u B ) + N S AB (u A T a lfl t A )(i B T a ^u B )] 

AytB 

where T a = A a /2, {A,B} = {L,R}, and the exponent 1 and 8 denotes a color singlet and a color octet 
interaction, respectively. 



III. RESULTS 



A. Parameter sampling method 



In order to find the best set of parameters that fits the data we use a Markov Chain Monte Carlo 
(MCMC) approach. We start with a random point from the multi-dimensional parameter space for the 
chosen model. The x 2 f° r this point is calculated and a likelihood is assigned to it. This likelihood is a 
measure of how well a set of data is reproduced for a given point in the model parameter space. Following 
the notation of |4Jj, this function is defined as 



G(0,O exp ,AO) =exp 

where 



- x 2 (0,O exp ,AO) 



(13) 



and O is the value of the observable for a given point of parameter space, O eX p is the central value of 
the observable and AO is the la error. The absolute value of the likelihood function is irrelevant for 
our analysis. What is relevant is the ratio of likelihoods, that is, the comparison of the likelihoods of 
two consecutive points in the chain. We follow the Metropolis-Hastings (MH) algorithm - the Markov 
chain is started from a random initial value in parameter space with a given likelihood that depends on 
the constraints imposed by the data set. Next, a new point is generated randomly with a probability 
distribution centred around the old point and the likelihoods of the two points are compared - if the 
likelihood of the next point is larger than the one for the current point the next point is appended to the 
chain, otherwise the current point is replicated in the chain. We have repeated this procedure with 10 
different random starting points. We scan with flat priors (i.e. a linear sampling over the parameters) 
and we have checked that the chains have a very good convergence behaviour. In all calculations of the 
top production cross sections we use the top mass as the renormalization and factorization scale. Wetake 
m t = 175 GeV and to take into account the NLO corrections we have chosen a A:— factor of 1.41 
Further, we use a bin-wise scaling for the m t i distribution to emulate the m t i dependent k— factor. 
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TABLE III: Relation between the constants presented in the plots and the independent FCNC operators in the 
strong and electroweak sectors. 



B. Strong and Electroweak operators 



In order to simplify the notation we have replaced the original constants from the operators in the strong 
and electroweak sectors by cvj with i = 1, 8. This correspondence between constants and the operators 
themselves is presented in table Hill As an example, the first operator C^<4 wou ld appear in the effective 
Lagrangian as 



x uG<p 

A 2 



/r\ut 



(15) 



Considering A = 1 TeV then a\ is defined as 



a. 



uG <t> TeV 2 



A 2 



(16) 



which renders a\ dimensionless. Similar definitions hold for the remaining on constants. Table flTTI shows 
the relations between all the constants shown in the plots and the independent FCNC operators in the 
strong and electroweak sectors. 




FIG. 5: The x 2 as a function of ai, 03 and o?7 with each parameter taken non-zero at a time. The most favoured 
values are ai = and |ct3 1 7^ 0. There are no preferred values for 07. 



We first present our results for the Strong and Electroweak FCNC operators (SEFCNC). In Fig. [5] we 
present the \ 2 as a function of at, a^ and a.7, keeping only one of the coefficients non-zero at a time. 
These three curves are representative of the x 2 distribution behaviour for the complete set of SEFCNC 
operators. In fact, we can group operators a\ and ai as for both a\ — and ai = are the most favoured 
values. The operators that are preferentially non zero when taken one at a time are a^, CI4, a§ and a^. In 
this case the preferred values are close to ai — ±4 (see a^ in Fig [5]) . Finally both a^ and as seem to be 
completely unconstrained as they have an almost flat x 2 distribution for the entire ai range presented. We 
have then proceeded to scan over the 8 parameters (pn, i = 1 — 8) using the MCMC method with flat prior 
as described in the previous section. The range for all parameters was chosen to be — 10 < ai < 10. The 
complete set of 14 experimental observables, presented in the introduction, is used to calculate the x 2 an d 
hence the likelihood. After the likelihood mapping for the model, we have obtained the one dimensional 



FIG. 6: One dimensional likelihood distribution of the parameters ai to as after the fit. 



likelihood distribution of the parameters which is presented in Fig |6] It is clear from the figure that both 
ai and «2, the FCNC operators stemming from the strong sector, are strongly constrained to be in the 
range —2 to 2. Operators 03 to ag, the chirality- flipping FCNC operators coming from the electroweak 
sector, have to be in the range —4 to 4. Finally the chirality-conserving operators from the electroweak 
sector, 07 and as are very mildly constrained and, as we will show later, the center-peaked shape of 
the distribution is only a reflection of the correlations of these parameters with the constrained ones. In 




6.5 7 7.5 0.1 0.15 0.2 0.2, 0.3 0.4 




0,02 0.04 0.06 0.06 0.08 0.1 0.120.140.16 0.1 0.2 0.3 

A FB( m ,r> A F B ( Y ,< 1 > A F B ( Y ,> 1 > 

FIG. 7: One dimensional likelihood distribution of the total cross section and all asymmetries after the fit. 

order to understand how well the model fits into the experimental observables we present in fig [7] the one 
dimensional likelihood distribution of the observables after the fit. Considering the central values of the 
observables presented in the introduction, we conclude that the SEFCNC set of operators prefer a lower 
value of the cross section while generating suitable values for the total forward-backward asymmetry. The 
values for Apg^m^) and A t FB (y^') are below their experimental central values but well within the error 
bands. This shows that there is some compromise between the values of the parameters in the attempt 
to fit all observable simultaneously giving rise to a slight difference between the input observables and 
the ones originated from the posterior probability distribution functions (pdfs). 

We now move to the study of the possible two-dimensional correlation between pairs of parameters. 
In Fig. [5] we present the two-dimensional correlation plots for the most representative scenarios. It is 
clear from the figure that there is no correlation between ct\ and 02- Furthermore, these operators are 
very strongly constrained. On the other hand, there are several pairs of values that cannot be zero 
simultaneously. This is the case of (03,04) - the ones from SU(2), (a 5 ,a 6 ) - the ones from U(l) and 
(03,05), (04,05) - these are the U(l) and SU(2) combination where the indices of the operators O l i are 
the same as for example in 0^ w ^ and O^^. For the pairs (04,05) and (03,05) the preferred values 
lie in the region o| = a\ and o| = a\ respectively. This happens to the combination of SU(2) and 
U(l) operators with the ij indices exchanged. Finally operators 07 and as do not appear to be much 
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FIG. 8: Two-dimensional correlation plots for the strong and electroweak FCNC parameters. Only the most 
relevant correlations are shown. The shaded areas are the ones where the values of on reach their highest 
probability - the darker ones represent 95 % CL while the lighter ones are for 68 % CL 

constrained when plotted against the remaining operators. There are however mild correlations - if we 
take for instance the pair (a4,ay) it is clear that for < 0, ay prefers to be positive and if ay > 0, ol-j 
prefers to be negative. 




FIG. 9: Likelihood distributions for the most relevant combination of parameters. 

With the hints from Fig. [5] about which parameters prefer to be non-zero after the fit, we have tried 
to understand if one could make a more strong statement about the appearance of new physics related 
to the Strong and Electroweak dimension six FCNC operators. We note that the contributions of 0:7 
and as are irrelevant because the change in likelihood is very small when these parameters are varied as 
shown in figures O and [51 On the other hand, a\ and a2 can lead to a large change in the likelihood - 
the preferred points are therefore ot\ = = ct2- Hence, we look at the most relevant combinations of the 
remaining parameters. The likelihood distributions for those combinations are shown in Fig. [9] It is clear 
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that all the correlated pairs of parameters prefer to be non-zero simultaneously, like for instance (03,04), 
which have a peak between 5 and 10. Again, the likelihood plot for a 2 — a| peaks at indicating that 
a| = a§ is the preferred parameter choice as also seen in Fig. El However, the most interesting case is 
the likelihood for a 2 + a 2 + a 2 + ctg - in this case we are certain that at least one of the four parameters 
has to be non-zero in order to fit the data. This is a very strong statement because it means that new 
physics coming from these operators can help curing the asymmetry discrepancy and in order to solve it 
at least one of the operators has to be present. 
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FIG. 10: Two-dimensional correlations between cross sections and asymmetries and between the different asym- 
metries. 



In Fig. [TU]we present the two-dimensional correlation between several observables after the fit. In the 
first row one can see that there is a negative correlation between asymmetries and total cross section. 
Hence, to get the right asymmetries the cross section moves to its lower preferred value. On the other 
hand, all asymmetries have positive correlations and are highly correlated - if one of them increases the 
other increases as well. Therefore, there is a tension between cross sections and asymmetries that reflects 
the difficulty of fitting all the observables with the set of SEFCNC operators. Nevertheless, a non zero 
contribution from the operators to ae provides a better fit than the SM one. 

In Table. IIV1 we show the best fit values along with 68% and 95% Bayesian confidence intervals 
(BCI) for all the parameters and selected observables. The BCIs are derived from the one-dimensional 
marginalized distributions, as shown in figures [5] and [6l while the best fit point is the one with least 
X 2 = 14.2. Thus, the best fit point does not need to be at the center of the marginalized BCIs. For 
the SM we have x 2 — 24.0 and it is the sizeable contributions from as, a§ operators that lead to the 
reduction in the x 2 for our fits. We again note that the combination a 2 + a 2 + a 2 + a\ > 7.5 with 97.5% 
CL, i.e. it is almost certainly non-zero. 

Wc have also listed the posterior BCI for the cross section and the asymmetries in table HVl The best 
fit value of the total cross section, and also the 95% BCI, are somewhat smaller than the measured 
central value. The same trend is observed for all the asymmetries except for the integrated asymmetry 
A l p B which is correctly reproduced and the A t FB (m t t — ) asymmetry which is most likely positive in our 
model. As previously discussed, the reduction of the cross section values and asymmetries is a result of 
the negative correlations between them. 
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18.91 
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TABLE IV: Best fit values and the Bayesian confidence intervals (BCI) for parameters and the observables. 




FIG. 11: Likelihood distribution for the parameters of the four-fermion Lagrangian, after the fit. 

C. Four fermion operators 

We now turn our attention to the four-fermion (4F) effective Lagrangian. We should start by mentioning 
that recently |42j, a complete calculation of the forward-backward asymmetry and of the total cross 
section of top quark pair production induced by 4F-operators was performed for the Tevatron up to 
0(a 2 s /h?). The results show that next-to-leading order QCD corrections can change both the asymmetry 
and the total cross section by about 10%. As discussed in section Hi CI there are a total of 12 independent 
operators for the study of ti production and under the conditions described previously which mainly 
means we are only considering the u-quark contribution in the initial state. We have scanned linearly 
over the 12 parameters from the 4F-Lagrangian using the MCMC method. The range chosen for all 
parameters was again from —10 to 10. In Fig. [11] we present the likelihood distribution for all the 4F 
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parameters, after the fit. A few comments are in order. First, operators in one row can only interfere with 
parameters in the same row. Second, only parameters in the first row interfere with the SM Lagrangian 
and consequently the main contribution for the asymmetry has to come from the parameters presented in 
the first row. This is clear from the plot as the four distributions in the first row are the only asymmetric 
ones - all other parameters in the following two rows have not only symmetric distributions but they 
show that the preferred value of these parameters is zero. However, in the case of 4F operators the 




2 4 6 8 2 4 6 8 -4 -2 2 0.2 0.4 0.6 

c c c c 

12 22 3 4 



FIG. 12: Likelihood distribution for the total cross section and for the asymmerties and for the six independent 
combinations of the 4F Lagrangian parameters, after the fit. 

cross sections and the asymmetries depend only on six combinations of the parameters. Therefore we 
have decided to present in Fig. fTJ] the likelihood distribution for those combinations together with the 
likelihood distributions of the total cross section and a few selected asymmetries. The relation between 
the new parameters and the original ones present in the 4F Lagrangian is 



Ci = 








c 2 = 


r 8 + r 8 

^LR ' ^RL 






Cl2 = 


- (cD 2 + (c RR ) 2 + 


9 
2 


[(CD 2 + (c RR ) 2 ] 


C22 = 


- (cIr) 2 + (c RL ) 2 + 


9 
2 


[(CD 2 + (CD 2 ] 


c 3 = 


°LL°Lfl T ^RR^RL 


+ 


2 [ClDlr + C RR C RL 


Ci = 


(Nl R ) 2 + (N RL f + 


2 
9 


[(nD) 2 + (Kl) 2 ] ■ 



It is clear that the experimental observables are well described by the fit. Regarding the parameters, the 
most relevant fact, that could already be inferred form Fig. frjl is that C 2 prefers to be non-zero and, for 
the same reason, the likelihood of both C\2 and C22 peaks at 1. A similar trend can now be seen in the 
two-dimensional correlations presented in Fig. [T21 It is clear that at 95 % CL the value zero is excluded 
in the top right plot. In the top left plot the value zero is still inside the 95 % CL contour. Regarding the 
correlations between cross section and asymmetries, and between pair of asymmetries, after the fit, the 
general trend is very similar to the one presented in the previous section for the strong and electroweak 
FCNC operators. Therefore we will make no further comments on those correlations. 

In Fig. [M] we present two dimensional correlations between C\ and C2 and the total cross section, total 
asymmetry and the parameter 622- We see that while C\ is positively correlated with both the cross 
section and the asymmetry, C2 is positively correlated with the cross section but negatively correlated 
with the asymmetry. Furthermore C±, C2 and C22 are all negatively correlated with each other. Finally, 
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FIG. 13: Two-dimensional correlations between the parameters that can give a significant contribution to the 
asymmetry. Also shown are typical examples of the correlations between cross section and asymmetries and 
between two asymmetry observables. 




6 1 1 I 1 -3.5 I J 

-10 12 -10 12 -10 12 

c c c 

1 1 1 




-3 -2 -1 -3 -2 -1 -3 -2 -1 

c c c 

2 2 2 



FIG. 14: Two-dimensional correlations between the parameters Ci and Ci and the total cross section, total 
asymmetry and the parameter C22. 



we conclude that either C2, C22 or both have to be non-zero which is not surprising given the relations 
presented in eqs. (3.5). 

In Table. [Vj we show the best fit point along with 68% and 95% BCIs. The best fit point is the one 
with least \ 2 — 6.28. As seen in figure fT3l only C\ B operators have relevant contributions to both the 
asymmetries and the cross sections. The weak operators, C\ B do not interfere with the SM diagrams, 
contributing therefore more to the cross sections and much less to the asymmetries. Thus, they are 
strongly constrained through the measured values of the cross sections. The N l AB operators contribute 
only to the cross sections and consequently are also strongly constrained and irrelevant as possible new 
physics contributions. 

Again, due to the negative correlations between the cross section and the asymmetries, there is a 
slight tension in the fits. This leads to a mild preference for lower values of the total cross section. The 
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A FB (Y t < 1) 


0.121 
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0.128 
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A FB (Y t > 1) 


0.420 


0.281 
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TABLE V: The table of best fit values for the 4F case along with 68% and 95% BCI. 



asymmetries, on the other hand, are reasonably well reproduced. We note that, ^4^ B (rn tt — ) prefers to 
be positive with 4F operators. 



IV. BOUNDS ON THE EFFECTIVE OPERATORS 



In this section we discuss all possible bounds on the dimension six effective operators described in the 
previous sections. Our goal is to ascertain whether the values of the couplings multiplying each operator, 
that could explain the measured asymmetry discrepancy, are still allowed by the available experimental 
data. We start by considering the dimension six FCNC operators. In section III B we have concluded that 
only effective operators stemming from the electroweak sector were likely to fit the Tevatron data on the 
top quark better than the SM. In fact, it is the sizeable contributions from the 03, ...,a6 operators that 
lead to the reduction in the \ 2 f° r our fits. We again note that the combination a§+a4+ag+a| > 7.5 with 
97.5% CL, it is almost certainly non-zero. Therefore, we have now to focus on the bounds for operators 
a3 to as to understand if such a high value of the constants is not in contradiction with experimental 
data from other sources. 

A very complete analysis on the electroweak FCNC operators was performed in [43| using not only all 
available data from B physics but also the data from direct FCNC top decays (the later will be updated 



in this work) . The bounds obtained on the operators taken one at a time are [43[ a§ < 0.81, a\ < 0.011 
and Oq < 0.096 while the best bound on a$ was shown to come from the direct constraint on BR(t — > qZ) 
and BR(t — > 57). Therefore, to satisfy a\ + a\ + a\ + a\ > 7.5 one needs ctg ~ 6.5. However, such a value 
of 0:5 would imply that BR(t — > qZ) « 3.7% and BR(t — > qj) « 6.3%. The most recent direct bounds on 
BR(t qj) and BR(t qZ) arc the ones from the Tevatron, 3.2% [Hi and from the LHC, 1.1% [H, 
respectively. Hence, it is clear that such high values of as are disallowed by Tevatron and LHC data on 
the direct searches for FCNC top decays with a photon or a Z-boson in the final state. Furthermore, 



Other analysis based on B physics observables and electroweak precision constraints were also performed in [3 leading 
to similar conclusions. 
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indirect bounds from HERA, where bounds on cross sections are converted on bounds on the branching 
ratio, set a limit BR(t — > (77) < 0.5% [13] ■ Also, a combine study on B physics and Tevatron data on 
top quark production cross section places an indirect bound on the sum of the FCNC branching ratios 
forcing them to be below the percent level [48| . In conclusion, experimental data from very different 
sources constrain the operators that could explain the asymmetry in such a way that we consider that it 
is very unlikely that the observed discrepancy could be explained by these operators. 

Contrary to the the dimension six FCNC operators, there are no useful bounds on the four fermion 
operators involving two top quarks and this is even more so if the top is right-handed. Therefore, only 
the LHC could place constraints on these operators. However the values of the constants C, and Cij 
that could help explain the discrepancy give an extra cross section that is always below 10 pb even for 
y/s = 14 TeV. Hence, given the error of tt production cross section it is very unlikely that these operators 
will be constrained in the near future. 



V. DISCUSSION AND CONCLUSIONS 



In this work we have used a dimension six Lagrangian with FCNC interaction together with four-fermion 
operators to gain some insight in understanding the discrepancy between the experimental values obtained 
for the top pair production asymmetry and the corresponding SM predictions. We have build a minimal 
set of operators and we have used an MCMC approach to find the best simultaneous fit of all independent 
operators to the available data. Our conclusions regarding which operators give the best fit are as follows 

• Strong FCNC operators with coefficients 0:1,2 are most likely close to zero; 

• regarding Electroweak FCNC operators with coefficients 03 to a§ we conclude that one of them 
must be non-zero; 

• Electroweak FCNC operators with coefficients 017,8 are not relevant; 

• bounds on electroweak FCNC operators reveal that it is very unlikely that the asymmetry can come 
from new physics described by these operators; 

• Four-fermion operators with coefficients C\ B contribute to the asymmetries as the ones with co- 
efficients C\ B give small contributions; the 4F combinations Ci, C2, C12, C22 contribute to the 
asymmetries; 

• Four-fermion operators with coefficients N\ B contribute to the cross sections only; 

• there is in all cases some tension between cross section and asymmetries when a simultaneous fit 
to all observables is performed; 

• the contribution of 4F operators to the cross-section at LHC7 is of the order ±1.5 pb, which is 
allowed by the present estimates of the cross-section [49j . 

It is important to ask how do we figure out which operators are actually responsible for the asymmetry. 
To that end we note that the asymmetry, although called forward-backward, is actually a C-odd and 
for CP conserving interaction that can therefore be created by -P-odd interactions as well. Further, 
our operators contribute to the asymmetry in two ways: kinematically and dynamically. The i-channel 
diagrams with FCNC interaction naturally originate more top quarks in the direction of the incoming 
it-quark leading to a positive asymmetry as measured by the CDF collaboration. This coupling does 
not need to be chiral to produce the required asymmetry, although our operators are chiral. For the 4F 
case, there is no such kinematical asymmetry generation and it is dominantly generated by the unbalance 
between left and right chiral operators interfering with the SM contribution. Thus, in both cases, we 
have the presence of chiral interaction, which also incarnates in the form of polarization of the produced 
top-quark. Hence, a study of such polarization effects (Hoj as a function of rapidity will be able to 
provide a probe of possible new physics. Further, our operators are also constrained by the B-physics 
observables and for simplicity we have not accounted for them in our MCMC. We have nevertheless used 
the constraints from B-physics to conclude that FCNC operators are unlikely to account for the measured 
asymmetry at the Tevatron 
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To conclude, we remark that the A" B observed at CDF can be casted in terms of dimension-six oper- 
ators and we need more observables, from top-polarization and B-physics, to constrain them due to the 
multitude of these operators. 
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